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We derive the formal expressions needed to discuss the change of the twist-two parton distribution 
functions when a hadron is placed in a medium with relativistic scalar and vector mean fields. 



I. INTRODUCTION 

Presently, the quark distributions of the free proton 
are quite well known over a wide range of x and Q 2 . 
However, the theoretical understanding of these distri- 
butions is somewhat limited. Although we know how to 
extract the quark distributions from the measured struc- 
ture functions within a NLO QCD analysis, we do not 
have sufficient control of non-perturbative QCD to cal- 
culate these distributions from first principles. Neverthe- 
less, we do know how to formulate the problem in terms 
of the proton matrix elements of certain local operators 
(see, for instance, Ref. |l(), which correspond to the mo- 
ments of the measured parton distributions. There has 
been considerable progress in calculating at least the first 
few moments using lattice QCD 2j - albeit at relatively 
large quark masses . 

When it comes to nuclear structure functions one needs 
to evaluate the matrix elements of these same operators 
in the nuclear ground state - a priori a much more dif- 
ficult problem. On the other hand, one knows that nu- 
clear structure functions (divided by the number of con- 
stituent nucleons) lie within 10-20% of the free nucleon 
structure function (except near the kinematic boundary 
for the free nucleon) Q|. It therefore seems reasonable 
to tackle the problem by computing the corrections to 
the structure function of a bound nucleon and then al- 
lowing for Fermi motion. In particular, one could con- 
sider as a starting point the case of infinite quark matter. 
Even this presents serious theoretical challenges, because 
one knows from numerous studies that one encounters 
large scalar and vector mean-fields in this problem 
and there has been no discussion of the formal aspects of 
parton distributions in such an environment (analogous 
to the discussion of Jaffe I] in the free case). 

We therefore begin with a formal development of the 
parton model for the case of a "proton" embedded in 
constant scalar and vector fields. Since we are interested 
in modelling QCD, asymptotic freedom will be imposed 
by hand, in that between the two hard collisions which 
define the forward Compton amplitude (for the leading- 
twist parton distributions) the quark struck by the pho- 
ton will be treated as free. The result of this formal 
investigation is of general interest as the formal proper- 
ties of the parton distributions, in terms of support and 
reflectivity or crossing symmetry (x — > —x), are guaran- 



teed. Even so, it is possible to make suitable definitions of 
both the valence and sea-quark distributions. In section 
II we derive the equations for these quark distributions 
in a bound proton. In section III we verify that they are 
normalized. Section IV is devoted to the investigation of 
the effects that the nuclear medium exerts on the distri- 
butions, while in section V we apply our results to the 
case of infinite, isospin symmetric, quark matter. Finally, 
in section VI we suggest directions for future work. 

II. QUARK AND ANTIQUARK 
DISTRIBUTIONS 

Our aim is to write down the formal expression for 
the quark and antiquark distributions of a bound proton. 
The in-medium proton momentum is denoted by P* — 
(P*°, P*), where the star superindex means that nuclear 
interactions have produced an effective mass and energy 
for the proton. As in deep inelastic scattering there are 
two independent variables to build the hadronic tensor. 
We will use P* and q, the photon momentum probing the 
nuclei, as those variables. Hence, the hadronic tensor for 
the bound proton is written as: 

+ ( P * ?L± Q ) ( p* ?-i q j jgVig!) 

(1) 

where Ff p (x*,q 2 ) and F BP (x* ,q 2 ) are the structure 
functions for the bound proton, and x* = —q 2 /2P* ■ q 
is the fraction of the bound proton momentum carried 
by the quarks. In the parton model, the structure func- 
tions are written as: 

2x*F 1 BP (x*,Q 2 ) = F BP (x*,Q 2 ) 
= E eUqr^l + qr^D + OiasiQ 2 ))^) 

i—u,d,s,... 

where Q 2 — ~q 2 , qf p (x*)(qf P {x*)) is the quark (anti- 
quark) distribution in the bound proton and 0(a s (Q 2 )) 
are the QCD corrections to the parton model. 
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In its simplest form, the quark distributions in the par- 
ton model are calculated from the handbag diagram. In 
the light-cone gauge, A + = 0, the distributions can be 
written as 



p*+ 

,BP ( ^ _ f , dz - 



*P*+z 

in i 

An 

< p*\i>Uz-)Mo)\P* >c |,w=o, (3) 



where P* + = (P*° + P* 3 )/y/2 is the plus component 
of the bound proton momentum and \P* > is the state 
vector of the bound proton. If the quark field operators 
in Eq. are expanded in terms of free plane waves, 
q BP (x) can be rewritten as: 



2tt J (2tt) 3 
<P*\bf(t)bl(0)\P* > c | zW=0 
P*+ f d 3 k 



-V a p(k) / dz 



ix* P*+ z -\-ik-z 



2tt J (2tt) : 
< P*|d?(t)df (0)|P* > c |, w=0 , 



(4) 



where 



U aP {k) = u\ a) {k)[l + 7°7 3 H/3)(fc), V aP {k) 



rization in- 



v (a)fo)[l + 7 7 ] v (P)(k), with a, (3 the pola: 
dices and the Einstein convention for repeated indices is 
understood. The time dependence of the creation and 
annihilation operators is calculated from: 

with H the QCD Hamiltonian operator. We assume that 
the state vector \P* > is an eigenstate of H with eigen- 
value P*°. Hence, the insertion of a complete set of in- 
termediate states in Eq. (0J implies that: 

q BP (x*) = 



2 

L 

n 

p*+ r d 3 k 



^ < P*\b a J\n >< n^JP* > c \ z+=z ^ 



2?r J (2vr 



■V a p{k) I dz e 



< P*W\n >< n\d P J\P* > 



■ ■'" ; -i(P*°-P*°)t+ik-z 

(6) 



where P*° is the eigenvalue of the intermediate state \n > 
after the action of the Hamiltonian operator. The inte- 
grals in z~ in Eq. © can be done: 



, lx 'P' + z -+i{P' a -P^)t-ik-z\ 

az e \ z +=z ± =0 

- 6(x*-k*+/P*+), 



p*+ 



(7) 



2tt 



c* P* + z- -i(P*° -P*°)t+ik-z\ 



S(x* + kl+/P*+), 



(8) 



where k* + = (P*° - P*° + k 3 )/y/2 is the plus component 
for quarks and k^ + = (P*° — P*° + / \[2 for antiquarks. 



We then have: 
q BP (x*) 



d 3 k 



U a0 (k)6(x* -k*+/P*+) 



(2tt) 3 

J2 < p *\bf\n >< n\b^\P* > c 

n 

rl 3 k 

— V aP {k)5(x*+k^/P* + ) 
J2<P*\d1\n><n\df\P*> c . (9) 

n 

We note that the second term contributes to the quark 
distribution because an antiquark with a negative mo- 
mentum fraction can be interpreted as a quark with pos- 
itive momentum, fc* + = —kl + . 

If we assume that quarks in the nuclear environment 
feel scalar and vector mean fields then we can write the 
vector potential as a shift in the quark and antiquark 
energies: 



(10) 



where V + is the plus component of the vector potential, 
and (k^) the plus component of quark (antiquark) 
momentum with masses modified by the scalar potential. 
Hence, from delta functions in Eq. © we have: 



x*P 



*+ 



xP 4 



ill) 



where x = fc+ /P + is the fraction of momentum carried 
by a quark inside a proton immersed in a medium with 
scalar but with no vector mean field. With the help of 
Eq. (|llfl we can isolate the effect of the vector potential 
on the quark distributions: 



q BP (x*) = (P*+/P+)q BP (x), 



(12) 



with the bound quark distribution without the vector 
potential defined as: 



HP 



(x) = — dz 



cP+j 



Pj 

-In 

< P*\^Uz-)^ + {Q)\P* 



=0.(13) 



As the effect of the vector mean field only changes the 
phase of the eigenstates, we conserved the notation \P* > 
for the bound proton with no vector potential when defin- 
ing q BP (x). Finally, because the scalar potential can be 
absorbed in the proton mass, it follows that the whole 
analysis of Jaffe [l| for the support of the Bjorken xgj 
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can be translated to the present case, meaning that the 
support for x is < x < 1. Similar relations to Eqs. Hll|) 
and (|12|l were also obtained by Mineo et al. @ , where the 
vector potential was treated as a gauge transformation. 

On the other hand, the antiquark distribution in the 
parton model is given by the following expression Q: 



q BP (x*) 



p*+ 

ill 



dz e 



-ix* P* + z 



< P*\^ + (0)i, + (z-)\P* > c | z+=z x =0 . 

(14) 

Inserting a complete set of intermediate states and per- 
forming the integrals over the z~ variable, as in Eqs. (JJJ 
and JSJ , we find the following expression for the antiquark 
distributions: 



-BP 



d 3 k 



U af3 (k)6(x* + k*+/P*+) 
Y J <P*T%\n><n\b^\P* > c 

n 

d 3 k 



:V af3 (k)5(x* 



kl + /P*+) 



(2tt) 3 

J2<P*\d^\n><n\d1\P* > 



(15) 



The expressions for the quark and antiquark distribu- 
tions, as given in Eqs. © and (fT5|l . have (as they should) 
the crossing symmetry: q BP (x*) 



-q BP {-x*). 



III. NORMALIZATION 



The quark distributions of the bound proton must 
have the correct normalization. That is, the integral of 
q BP (x*) — q BP (x*) over the allowed x* range must give 
the number of valence quarks of the bound proton. To 
this end, we will use Eq. (|llfl . together with the allowed 
support for x. It follows that the maximum value for x* 
is 1 — 2V + . while the minimum value of a;* is V + : 



l-2V i 



v+ 
p*+ 

~2^~ 



[q BP (x*)-q^(x*)]dx* 



dx*dz e 



-IX* P*^ 2 



=<P*|Vi(0hMQ)|P* > c , 



<p*hft(OiMo)|P* >c 

(16) 



which gives the quark number in the state \P* >: 

rl-2V+ 



[Q 



BP i 



-r.BPl 



V+ 



IV. 



q" (x*)]dx* = N BP -N± 



BP 

q " 



(17) 



THE IN-MEDIUM EFFECTS ON THE 
DISTRIBUTIONS 



The quark and antiquark operators appearing in the 
distributions @ and i|15fl are those of free quantum fields. 



The quarks which build the proton state are, however, 
not in free space: they are confined in a proton state, 
while the proton state is itself immersed in a nuclear 
medium. Our state vector has to be built from these 
bound quark operators. We shall denote by bp (dt) the 
annihilation operator of a quark (antiquark) in the bound 
proton. Let the bound proton state be written as: 



1^ >-F[ q ;l, q ;l q ^w >, 



(18) 



where P[g^,gt^,g^|] is a functional of the bound quark 
operators, |0* > is the effective vacuum (as seen by 
quarks in a nuclear medium) where the quarks bound 
in the proton live, and pi,p2,P3 are the individual mo- 
menta of the three valence quarks. Although only the 
valence quarks appear explicitly in the functional, it is 
understood that it may be populated by quark - anti- 
quark pairs and gluons: the notation only shows the net 
number of quarks and antiquarks. 
The effective vacuum is defined by: 



fe*|0* >= , d*J0* >= 0. 



(19) 



With this definition, we can calculate the action of the 
free quark and antiquark operators on the bound proton 
state. To do this, note that the quark field operator in 
Eq. J2Jl could have been expanded in any basis 01 ■ Sup- 
pose that we know the solution of the Dirac equation 
for the interacting theory describing the bound proton 
state, with solutions ubp(p,x) {vbp(p, x)) for the posi- 
tive (negative) energy part. The field operator expanded 
in this basis is written as: 

. . (20) 

with u'u = 1 [5J. If we compare it with the expansion in 
terms of the free fields, we get: 



d 3 p 
(2^)3 



dp i0 (t)B^(k,p)], (21) 



where: 



A a0 (k,p) = 



d 3 xu {a) '(k) l0 u^ ] P {p,x)e 



ik-x 



(22) 



B a/3 (k,p)= / d 3 xu (Q) (/fc) 7 o<j,(p,x)e- lfc ^. (23) 



The calculation of the anticomutator between the free 
and bound operators gives: 

{bm,bf(t')} t=t ,=A«P(k,p), (24) 
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and similar for the other operators. Thus the action of 
the free quark annihilation operator in the proton state 
results in: 



bl\P*> = A^(k > p i )G[bp.,bp i ]'y\0*> 



(25) 



where G\b*J , b*JP is some function of the effective quark 

Vj Pi 1 

operators after 62 acts on F, and the index 7 indicates 
that the resulting two quark states are in coloured states. 
Integration over pi, pj and p\ is understood. Thus, the 
action of the free operator on the effective vacuum yields: 



6S|0* >= 



d 3 p \ - 

(2tt) 3 ^ 



B^{k,p)\q},, > ■ (26) 



From Eqs. Ij25(l and Ij26(l we see that there will be contri- 
butions from two terms in the quark distribution, Eq. 
one when we have two quarks in the intermediate state 
and the other when we have three quarks and one anti- 
quark in the intermediate state. For the case n = 2, we 
define the coloured state: 



\n = 2>=G[b]?,bp k y\0*>, 
and for the case n = 4: 

|n = 4>=F[6g,&g,6g]|^, 7 > 



(27) 



(28) 



In a completely analogous way to what was done in 
Eqs. (|21|l - Ij23|l for the free quark operator, we can ex- 
press the free antiquark operator in terms of the bound 
quark and antiquark operators. We will call C a °(k,p) 
the overlap analogous to the B a P(k,p) term, with the in- 
terchange between the u and v Dirac spinors. Similarly, 
D a P(k,p) will be the analogue of the A a ^(k,p) overlap, 
with the v Dirac spinor replacing the u Dirac spinor ev- 
erywhere. 

Using the states defined in Eqs. (|27() and i|28|) . the 
quark distribution Eq. is finally written as: 



BP 



(x*) 



d 3 k d 3 p 



U a p(k)6(x* -k* q +/P*+) 



7 



(2tt) 3 (2tt) 3 

(< n = 2\n = 2 > A ia \k,p)A lP (k,p) 



< n = 4|n = 4 > B' ra "'{k,p l )B^ l3 (k,p) 

|||f>^<£«z-+*r/p-+) 

"^2 (< n = 2 \ n = 2 > D' ia Hk,p)D'> i:l (k,p) 

7 

- <n = 4|n = 4>C 7at (fc,p')C 7/3 (fc,p)) 

(29) 

where integration over all the internal momenta of the 
quarks in the intermediate states is understood. 



V. THE QUARK MATTER CASE 

A central question from the point of view of nuclear 
physics involves the changes to the quark and antiquark 
distributions of a bound proton. Since one must develop 
a reliable model of both the free proton and the binding 
of nucleons starting from the quark level [8j , this problem 
is rather complicated. We intend to report on our inves- 
tigation of that problem in future work. For the present, 
we have chosen to illustrate the formal ideas developed 
here by applying them to a toy model, namely the quark 
distributions of isospin symmetric quark matter in which 
each quark feels a scalar potential, — V^, and a vector 
potential, V^. This is the premise of the Quark Meson 
Coupling (QMC) model @ which has been used success- 
fully to calculate the properties of nuclear matter as well 
as finite nuclei [Tol . Most recently it has also been 
used to derive an effective nuclear force which is very 
close to the widely used Skyrme III force (Except 
that in QMC the quarks are confined by the MIT bag, as 
well as feeling the mean-field scalar and vector potentials 
generated by the surrounding nucleons.) In the mean 
field approximation, the Dirac Equation for the quark in 
infinite quark matter is written as: 

[ij ■ d - (m - V?) - lvV«W QM {*) = 0- (30) 

Following Eq. (|20|l . we write the field operator in terms 
of the solutions of Eq. (|30|l : 



V'qm 0*0 



d 3 p 
(2^)3 



E 



b p, a UQ M (p)e 



-ipZ t-\-ip-x 



+ dpy Q t(p)e< t -^}, (31) 

where Ug"2f(p) an( i v< qm(p) are the in- medium Dirac 
spinors with mass m* = m — U/, energy E*, and where 
p*° = E* + V v q for quarks and pt° = E* - V v q for anti- 
quarks. 

We now calculate Eq. 112211 for the A factor, using k — 
{p%,k) for the free quark and p — (p*°,p) for the quark 
in infinite quark matter: 

A a \Kp) = 

J d 3 xu^(k) l0 u^ M (p)e^-^. (32) 

The integral over x can be done, giving a Dirac delta 
function, which implies that: 

A^(k,p) = {2*Y5{p- k)u^(k) loU ^ M (p). (33) 

Similarly, Eq. I|23|l for the B factor becomes: 



B^(k,p) = (2tt)^(p+ fe)uW(fe) 7Q ^(p). (34) 



5J<«)f 



,03) 



QM y 
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If we go back to Eq. I|21|) relating the free and bound 
quark operators, we find: 

+^")(fc) 70 ^(-fc)d: t - /3 (t). (35) 

A similar calculation for the antiquark operator gives: 



df{t) = v^\k) lQ u^ M {~k)b*^(t) 
+v^\k) l0 v i Q P l I (k)df(t). 



(36) 



An explicit calculation shows that D a ^(k,p) — A a ^(k,p) 

and C a °(k,p) = B a P(k,p), which implies that charge 
conjugation holds between the bound quark and anti- 
quark operators. 

Although the quark and antiquark operators in Eqs. 
(|33|) and (|3l)|l are time dependent, the expressions for the 
quark distributions involve products like b^b, which are 
time independent: 



blh = bUt)bM. 



(37) 



The bound quark distribution, calculated from Eq. 
is then: 



d 3 k 



U a/3 (k)S(x* - k*+/P*+) 



E q + mE*+m* 
2E q 2E* 



(E q + m)(E*+m*) 



<P*\b*J b% \P* > 

' k,a K,p' 



E q + m E* + m* 



x 04* • *Xa)(*J* ■ k X s) < P*\d*_Y^\P* > 



d 3 k 



(27T) 



3 V a0 {k)S(x* + k*+/P*+) 



E q +m E* + m* 
2E q 2E* 




< P*\d* r d*J IP* > 

1 k,a k,/3' 



x (x\* ■ ibca)(x& ■ kXs) < P*\b*J^\P* > 



VI. CONCLUDING REMARKS 



(38) 



We have developed the formal framework for evaluat- 
ing the change in the structure function of a hadron when 
it is bound in a nuclear medium. In particular, the for- 
malism is able to deal with the separate changes in the 
quark and antiquark distributions, while preserving the 
necessary sum rules. Since much of our information on 
the parton distribution functions of the nucleon does in 
fact come from nuclear data this is an especially impor- 
tant issue [TH . IT3 . PHI . In order to illustrate the formal- 
ism we considered the parton distribution functions for 
quark matter with mean field scalar and vector poten- 
tials. In the immediate future it will be important to 
include the effect of confinement, as for example in the 
QMC model dEil or even the NJL model [aE! with 
proper time regularization 0, 0] . 
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